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Abstract 


The  effect  of  finite  sample  thickness  on  material  elasticity  measurements  marlp 
using  an  Atomic  Force  Microscope  (AFM)  has  been  calculated.  The  model  includes  an 
elastic  layer  on  an  elastic  foundation  and  simulates  sample  indentation  under  an  applied 
load,  for  rigid  axisymmetric  tips  with  conical,  paraboloidal,  and  hyperboloidal  profiles. 
The  results  show  that  a  common  approach  to  estimating  elastic  modulus  from  force- 
displacement  curves  can  lead  to  significant  error  that  depends  on  the  units  of 
measurement.  A  method  to  estimate  this  error  unambiguously  and  correct  it  is  proposed. 
In  addition,  it  is  shown  that  elasticity  estunates  for  monolayer  thick  samples  using  the 
force-modulation  technique  have  substantial,  sample  thickness-dependent  error.  Local 
thickness  variation  can  result  in  misleading  contrast  in  force  modulation  images  for 
samples  that  are  several  nanometers  thick. 
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Atomic  Force  Microscopy  (AFM)  measurements  are  often  made  on  thin  film 
samples  with  thickness  ranging  from  few  nanometers  to  microns.  Quantitative  elasticity 
measurements  using  AFM  are  usually  interpreted  assuming  semi-infinite  sample 
geometry.  In  this  paper  we  model  the  effect  of  substrate  on  elasticity  measurements 
made  on  thin  films.  This  effect  is  particularly  important  when  Yoimg  modulus  mapping  is 
performed  on  samples  exhibiting  inhomogeneous  thickness  variation’^  and  when  a 
breakdown  of  macroscopic  models  of  elasticity  is  tested.  The  modulus  of  elasticity  for 
typical  polymeric  and  biological  samples  is  usually  several  orders  of  magnitude  smaller 
than  the  modulus  of  the  underlying  hard  substrate^  We  calculate  the  sample's  indentation 
by  considering  a  rigid  axisymmetric  tip  under  an  applied  load,  and  we  report  a  significant 
error  arising  from  the  common  assumption  of  infinite  sample  thickness.  We  consider 
conical,  paraboloidal,  and  hyperboloidal  tip  geometries. 

Presently  the  most  common  techniques  for  measuring  elasticity  with  high  spatial 
resolution  are  indentation  measurements  ^  and  force-modulation  experiments, 
which  are  performed  using  AFM.  These  measurements  are  interpreted  using  a  model  for 
a  paraboloidal  or  conical  tip  elastically  indenting  a  half-space  sample.  However,  sharp 
AFM  tips  often  have  a  conical  shape  with  a  roimded  apex  (manufacturers  usually  specify 
10  -  100  nm  radius  of  curvature).  Thus,  to  simulate  this  shape  we  choose  a  hyperboloidal 
tip  profile.  Figure  la  shows  the  shapes  of  tips  we  use.  Figure  lb  shows  the  geometry  of 
indentation  and  introduces  the  variables  and  parameters  we  use  later  in  this  article. 

The  elastic  modulus  for  a  semi-infinite  sample  can  be  estimated  using  the  load- 
indentation  dependencies  as  described  by  Sneddon. 
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Equations  U-c  «,„aspond  to  oonioal.  pantlxtloidal  and  hyperboloidai  shapes  of 
r^olution,  tnspe^ively.  Hete  F  is  load,  F-  indentation,  a-  contact  area  radius,  F  - 1^ 
radius  of  curvature,  o-tip  semi-vertical  angle.  F .  Young  modulus  <r-  Poisson  ratio, , 
fflustrated  in  Figum  1 .  Fmm  Equations  la-b  it  nnty  be  seen  that  a  linear  fit  in  log(F)  vs. 
log(^  coordinates  has  an  intercept  related  to  the  Young  modulus: 


(2a) 


^  (2b) 

We  first  note  that  using  Equation  2.  (conical  shape)  for  tips  with  rounded  apex, 

results  in  considerable  e^r.  For  example,  in  order  to  estimate  the  Young  modulus  with 

a  teor  of  two  (using  hyperboloidai  tip  with  20-  semi-verhcal  angle  and  10  nm  radius . 

curvature,  according  to  Equations  Ic)  the  maximum  indentation  should  be  forger  than  2 

pnr.  On  the  other  hand.  Equation  2b  gives  ^tor  of  two  error  for  indentations  that  are  c, 

400  nm  using  the  same  parameters.  This  is  anticipated  since  Equations  le  transfonn  int, 
la  when  Je«rron(<?),  and  into  lb  when  o^f!  cot(c4. 

It  IS  mtunively  clear  that  if  a  soft  sample  rests  on  a  hard  foundation,  then  the  tip 
starts  to  "feel"  substructure  when  the  radius  of  contact  area  is  compamble  to  the  sample's 
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thickness.  This  results  in  force-indentation  curves  that  are  non-linear  in  logarithmic 

coordinates.  Figure  2  is  an  expHcit  illustration  of  such  an  effect.  The  Figure  shows  the 

indentation  of  a  -30  nm  thick  spin-cast  sampie  of  coblock  (polystyrene- 

polyvmylpyridine)  polymer  on  a  microscope  coversiip  substrate.^  The  tip  radius  was 
approximately  60  mn. 

For  this  type  of  response,  Young  modulus  estimates  can  have  an  enor  that  varies 
with  the  units  of  measurements.  Indeed,  the  error  depends  upon  the  distance  between 
intercept  of  the  linear  fit  line  (#=/  unU)  and  the  line  eorresponding  to  the  conditions 
assumed  when  deriving  Eqtations  1.  As  illustrated  iu  Figure  2,  such  an  eiror  can  even 
change  sign.  Using  the  units  shown  on  Figure  2.  the  estimated  modulus  would  be  two 
times  smaller  than  the  modulus  obtained  when  a/h  value  is  small. 

To  clatily  this  problem  we  have  perfoimed  numericai  calculations,  which 

simulate  the  indentation  of  an  elastic  layer  bonded  to  an  elastic  foundation.  Our 

computations  are  based  on  an  elastic  layer  model  by  Dhaliwal  and  Ran  Force- 

indentation  cun^s  are  obtained  by  solving  a  Fredholm  integral  Equation  of  the  second 
kind; 


^  J  ‘  (3) 

where  K(:c.O  is  the  kernel  of  the  integral  Equation,  £  and  .rare  the  Young  modulus  and 
the  Poisson  ratio  of  the  layer,  respectively.  The  kernel  is  defined  as: 

£(A:,r)=  2  (2a)cos  ( ^lu]coJ S-x,,  U 

0  KhJ[h 
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n  is  the  ratio  of  Lame  constants  /i  =  (E/2(l  +  o-))/(^,/2(l+o-,)  in  Equation  (3) 

is  the  tip  shape  function,  and  here  we  use  the  following  shapes  : 


/(x)=axcota, 

fix) = *  cot  *  a|^i/(ax/ficotfl:y  +  l  - 1  j. 


conical  tip 
paraboloidal  tip 

hyperboloidal  tip 


(4) 


where  a  is  the  radius  of  the  tip-sample  contact  area,  R  -  radius  of  the  curvature  of  the  tip 
apex,  or  -  tip  semi-vertical  angle  as  shown  on  Figure  1. 

To  find  the  indentation  S,  we  use  the  condition  that  the  normal  component  of 
stress  remains  finite  around  the  circle  of  contact  between  the  tip  and  the  sample.  A 
solution  of  (3)  must  satisfy 

«>(O=0  (5) 

We  solve  Equation  (3)  numerically  to  determine  the  indentation.  For  each  value 
of  the  contact  radius  a  we  iteratively  find  the  indentation  5  such  that  the  solution  of  (3) 
satisfies  (5),  and  we  find  the  load  force  according  to  : 

0 

To  estimate  the  error  in  the  Young  modulus  unambiguously  we  plot  logarithmic 
force  indentation  curves  in  normalized  coordinate  units:  ** 


6 


s= — i — 

Acot(a) 

p-^g-o  ‘ 

2A*  cot(fl:) 
y  S-2R 

A" 


conical 

tip 


paraboloidal 

tip 


(7a) 

(7b) 


For  conical  and  paraboloidal  tips,  force  vr.  indentation  curves  plotted  in  these 
eoordinates  ehange  only  when  ratio  of  Lame  coefficients,  m,  changes,  n  shows  the 
mismatch  between  the  rigidity  modulus  of  the  top  layer  and  the  rigidity  modulus  of  the 
substrate.  When  /t=/  and  o*oj,  the  kernel  vanishes  and  Equation  3  becomes  equivalent 


to  the  Abel  integral  equation  considered  by  Sneddon  “.  The  elastic  layer  indentation 
curves  are  not  linear  in  logarithmic  coordinates.  Therefore,  the  slope  and  the  intercept  of 
straight  lines  fitted  through  the  data  points  depend  upon  the  maximum  indentation  value. 
The  Young  modulus  can  be  estimated  in  normalized  coordinates  using  Equations  la  and 

lb  with  force  and  indentation  substituted  firom  Equations  7a  and  7b  for  any  - - 

indentation  value.  Hereafter  we  use  the  term  Young  modulus  to  indicate  the  true  Young 
modulus  of  the  top  layer  material,  and  "estimated  Young  modulus"  to  mean  the  Young 
modulus  estimated  from  the  linear  fit.  We  calculate  the  ratfo  of  this  estimated  Young 
modulus  to  the  modulus  of  the  top  layer  (we  caU  this  ratio  a  reduced  Young  modulus)  and 
plot  It  versus  normalized  maximum  indentation.  Plots  for  different  values  of /t  are  shown 
m  Figure  3.  The  left  panel  shows  a  calculation  for  a  conical  tip,  and  the  tight  panel  shows 
the  calculation  for  a  paraboloidal  tip.  ^  values  are  shown  to  the  right  of  each  curve.  We 
used  equally  spaced  lvalues  in  our  calculatioa  In  the  AFM  experiment  this  is  collect 

when  the  cantUever  spring  constant  is  considerably  larger  than  the  sur&ce  spring 
constant. 
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The  reduced  modulus  from  Figure  3  can  be  used  as  a  correction  factor  to  obtain  a 
better  estimate  of  the  Young  modulus  for  soft  samples.  Indeed,  the  lines  for  //  less  than 
10-'  are  almost  indistinguishable.  Therefore  these  lines  can  be  used  as  correction  factors 
for  samples  with  Young  modulus  ~0.1  GPa  and  less  and  substrate  modulus  -100  GPa. 

For  larger  values  of  jj  an  iterative  procedure  can  be  used  to  correct  both  E  and  /r. 
As  an  explicit  example  using  the  data  plotted  on  Figure  2,  we  calculate  5  and 
P  accordmg  to  7b.  We  take  0=0.33  as  a  reasonable  suggestion.  The  estimated  Young 
modulus  IS  calculated  from  the  intercept  of  linear  fit  in  log-log  coordinates  according  to 
equation  E  =  3exp(A)/V2  ,  where  b  is  the  intercept.  Thus  we  obtain  £=133  MPa.  The 
maximum  normalized  indentation  is  =0.88  .  Using  Figure  3  we  can  conclude  that  £  is 
overestimated  by  a  fector  -2.1,  thus  the  true  modulus  £=63  MPa.  This  value  of  Young 
modulus  happens  to  be  close  to  the  one  used  to  calculate  the  dashed  line  on  Figure  2,  58 
MPa.  We  note  that  our  indentation  values  are  not  equally  spaced,  which  contributes  to 
inaccuracy  in  the  modulus  determination.  It  is  important  to  know  the  layer  thickness 
accurately  to  determine  the  correction  for  £  (since  normalized  indentation  in  7b  is 
proportional  to  h  ).  On  the  other  hand,  if  the  Young  modulus  is  known,  the  steep 
dependence  of  £  on  A  can  be  used  to  estimate  the  layer  thickness. 

Force  -  indentation  curves  for  a  hyperboloidal  tip  were  transformed  according  to 
Equation  7b.  The  results  are  shown  on  Figure  4,  which  shows  the  reduced  Young 
modulus  plotted  against  the  maximum  normalized  indentation,  for  /isO  and  different  h/(R 
cot  a)  ratios.  If  the  h/(R  cot  a)  ratio  is  small  (less  than  0.1)  then  the  hyperboloidal  tip 
mdentation  curves  are  very  close  to  the  paraboloidal  tip  curve  (dashed  line).  Also  small 
values  of  normalized  mdentation  correspond  to  a  small  ratio  of  maximum  contact  area 
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radius  to  layer  thickness.  In  this  case,  hyperboloidal  tip  indentation  featutes  are  similar  to 
paraholoidal  tip  indentation.  Odterwise,  hyperl»loidal  punch  indentation  ean  not  be 
described  by  simpie  models  and  should  be  treated  explicitly. 

In  addition,  we  wouid  like  to  comment  on  the  use  of  the  force  modulation  method 
.0  determine  tte  eiasticity  of  thin  layers.  In  this  method  the  sample  oseiibtes  at  a  cemin 
ftequency.  and  this  oscUlation  is  sensed  by  the  tip.  which  contaas  the  sample.  The 
amplitude  of  the  tip  osciliation  depends  upon  a  sampie  stifihess  and.  for  a  semi-inflnite 
sampie.  the  Young  moduius  can  he  estimated  using  the  foUowing  Equation 


E  = 


6RF 


(8) 


where  is  cantUever  spring  constant,  H  -  tip  radius  of  curvature,  f  -  average  ioad 
applied  to  the  sample,  z  -  modulation  amplitude  of  the  sample,  ,ip  oscillation 
amplttude.  For  a  thin  sample  such  an  estimate  ean  result  in  substantial  eiror.  Figure  5 
shows  the  reduced  Young  modulus  estimated  according  to  Equation  8  vs.  appUed  ioad  for 
tuples  of  different  thickness.  The  eaieulation  parameters  are  given  in  the  figum's 
caption  Figure  5  demonstmtes  that  a  change  in  the  sample  thickness  may  msuh  in  an 
apparent  softness  contrast  in  force-modulation  images. 

Figure  6  shows  the  same  dependence  as  Figure  5,  where  the  lines  are  calculated  in 
normalized  coordinates  using  Equations  7b  and  8  for  seveml  vaiues  of  ^  If  ,he  layer 
thickness,  tip  radius  and  elastic  parametem  of  sul»trate  o»  are  known  and  a 
reasonable  assumption  about  Poisson  ratio  of  elastic  layer  ean  be  nmde.  the  reduced 
elastic  modulus  from  Figure  6  can  be  used  to  cotr^n  the  error  in  measured  modulus. 
Since  we  do  not  know  correct  value  for  ^  (it  depends  on  true  Young  modulus),  an 
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iterative  procedure  can  be  employed  to  determine  the  correction  coefficient  E  .  Let  E*  be 
the  symbol  for  the  Young  modulus,  as  estimated  from  equation  8.  To  obtain  the 
correction  factor  for  E*,  the  corresponding  average  normalized  load  can  be  obtained  from 
Equation  7b.  This  gives  the  fixed  position  on  the  abscissa  of  Figure  6.  The  ratio  of  E*  to 
the  Young  modulus  of  substrate  (Ei)  gives  an  approximate  //  value.  From  Figure  6,  an 
approximate  correction  fector  E  can  be  fr)und,  from  which  we  frnd  a  better  estimate  for 
This  gives  an  adjusted  correction  fector  E,  which  in  turn  improves  the 
estimate  of  The  iterations  converge  in  a  few  steps  to  a  best  value  for  //,  resulting  in  the 
correction  fector  for  E*.  For  an  explicit  example,  we  consider  E*  =0.5  GPa,  £/=  100 
GPa,  and  F"  =1.  The  iterative  process  using  these  parameters  is  shown  with  circles  on 
Figure  6.  Initially  we  have  //«0.005.  From  the  figure  we  can  see  that  for  this  E  x3. 
This  decreases  the  estimate  of  ju  to  0.0016,  giving  us  E  »4.5  .  The  next  step  results  in 
^«5 ,  and  subsequent  iteration  does  not  change  E  significantly.  Thus,  we  can  conclude 
that  the  true  elastic  modulus  is  close  to  0.1  GPa. 

In  summary,  we  have  used  an  elastic  layer  on  an  elastic  foundation  model  to 
calculate  the  error  associated  with  using  elastic  half-space  model  to  estimate  the  Young 
modulus  in  AFM  indentation  and  force  modulation  experiments.  We  find  that  the 
application  of  semi-infinite  sample  models  can  result  in  unpredictable  and  significant 
error  in  Young  modulus  estimates.  Normalized  coordinates  should  be  used  to  obtain 
consistent  results.  Contrast  in  the  force-modulation  imaging  may  be  related  to 
inhomogeneity  in  the  sample  thickness.  We  have  proposed  a  method  for  reducing  the 
error  related  to  finite  sample  thickness. 
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Figure  Captions 


1.  (a)  AFM  tip  profiles  considered  in  this  paper,  a  -  semi-vertical  angle,  R  -  radius  of 

curvature  of  the  tip  apex,  (b)  Geometry  of  sample  indentation:  S-  indentation;  F- 

load;  a  -  contact  area  radius;  kc  cantilever  spring  constant;  h  -  sample  layer  thickness, 

E,  Ej,  a,  <Ti  -  Young  moduli  and  Poisson  ratios  of  the  sample  layer  and  substrate, 
respectively. 

2.  Deviation  of  load-indentation  from  power  law  dependence.  Force  -  indentation  data 
was  coUected  on  a  ~30  nm  thick  spin-cast  sample  of  coblock  (polystyrene- 
polyvmylpyridme)  polymer  on  a  microscope  coverslip  substrate.  The  tip  radius  was 
approximately  60  nm.  Fitting  to  a  straight  line  in  logarithmic  coordinates  gives  an 
error  m  the  Young  modulus  which  depend  upon  the  origin  that  is  chosen  for  log(^, 
which  IS  determined  by  the  units  of  measurement.  Here  a  solid  line  is  the  best  linear 
fit  to  the  force-indentation  values  that  are  shown  as '+'  symbols.  The  dashed  line  was 
calculated  according  Equation  la  with  the  Young  modulus  calculated  from  the  data 
subset  with  a  small  a/h  ratio.  The  brackets  indicate  a  mismatch  in  the  intercept,  which 
IS  related  to  the  error  m  the  Young  modulus  (Equations  2).  Two  possible  positions  of 
log(^  origin  are  shown. 

3.  X  axis  -  normalized  maximum  indentation,  Y  axis  -  ratio  of  Young  modulus 

estimated  from  Imear  fit  to  modulus  of  top  layer.  Left  panel  -  conical  tip,  right  panel 
-  paraboloidal  tip.  Different  lines  correspond  to  different  ratio  of  Lame  coefficients, 
lines a-g  //correspond to  1, 0.46,  0.22,  0.1,  0.046, 0.01,  and  0.001, respectively. 
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4.  Reduced  Young  modulus  ploKed  vs.  normalized  maximum  indentation  for 

hyperboloidal  tip  (soBd  lines).  Coordinates  tmnsformed  acconling  to  the  paraboloidal 

tip  Equations  7b.  Hem  .0^  and  lines  (from  the  top)  conespond  to  h/(R  co,  a) 

ratios:  1.7, 1.2,  0.84,  0.58, 0.36  and  0.18.  For  comparison  we  show  a  calculation  for 
the  paraboloidal  tip  (dashed  line). 

5.  Ratm  of  the  modulus  estimated  using  Equation  (8)  to  actual  value  given  as  function  of 
avemge  load  for  samples  of  different  thicicness.  Diffemnt  lines  correspomi  to  diffemnt 
layer  thickness.  The  thickness  (from  1  to  1 00  nm)  is  indicated  on  the  grapL  Other 

parametem:  tip  radius  of  curvature  10  nm,  yr-O.OI,  £=,  oPa,  cH).  k,  10  N/m,  sample 
oscillation  amplitude  0.02  nm. 

6.  Plots  of  the  reduced  modulus  vs  the  nonnalized  force  for  different  yr  values  illustrate 

a  possible  error  in  the  Young  modulus  as  demmnned  using  the  force-modUation 

technique.  Starting  from  the  top  yn  0.001,  0.0022,  0.0046,  0.01, 0.022, 0.046,  0.1, 

0.22  and  0.46.  Circles  show  an  iterative  process  to  determine  the  cormction  fector  for 
estunated  Young  modulus. 
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